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Abstract 

In hierarchical cosmologies the evolution of galaxy clustering depends both on cos- 
mological quantities such as £1, A and P(k), which determine how collapsed structures 
- dark matter halos - form and evolve, and on the physical processes - cooling, star 
formation, radiative and hydrodynamic feedback - which drive the formation of galax- 
ies within these merging halos. In this paper, we combine dissipationless cosmological 
N-body simulations and semi-analytic models of galaxy formation in order to study how 
these two aspects interact. We focus on the differences in clustering predicted for galax- 
ies of differing luminosity, colour, morphology and star formation rate, and on what 
these differences can teach us about the galaxy formation process. We show that a 
"dip" in the amplitude of galaxy correlations between z = and z = 1 can be an impor- 
tant diagnostic. Such a dip occurs in low-density CDM models, because structure forms 
early and dark matter halos of mass ~ 1O 12 M0, containing galaxies with luminosities 
~ L*, are unbiased tracers of the dark matter over this redshift range; their clustering 
amplitude then evolves similarly to that of the dark matter. At higher redshifts bright 
galaxies become strongly biased and the clustering amplitude increases again. In high 
density models, structure forms late and bias evolves much more rapidly. As a result, 
the clustering amplitude of L* galaxies remains constant from z = to z = 1. The 
strength of these effects is sensitive to sample selection. The dip becomes weaker for 
galaxies with lower star formation rates, redder colours, higher luminosities and ear- 
lier morphological types. We explain why this is the case and how it is related to the 
variation with redshift of the abundance and environment of the observed galaxies. We 
also show that the relative peculiar velocities of galaxies are biased low in our models, 
but that this effect is never very strong. Studies of clustering evolution as a function 
of galaxy properties should place strong constraints on models of galaxy formation and 
evolution. 
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1 Introduction 



Local galaxies are highly clustered. On large scales they are organized into a network of 
sheets and filaments which surround large underdense regions, usually referred to as voids. 
On smaller scales galaxies are found in gravitationally-bound groups and clusters. Accord- 
ing to the standard theoretical paradigm, the structures observed today were formed by the 
gravitational amplification of small perturbations in an initially gaussian dark matter density 
field. Small scale overdensities were the first to collapse, and the resulting objects subse- 
quently merged under the influence of gravity to form larger structures such as clusters and 
superclusters. Galaxies formed within dense halos of dark matter, where gas was able to reach 
high enough overdensities to cool, condense and form stars. 

In this hierarchical formation picture, the clustering of the dark matter, as measured by 
the amplitude of the matter correlation function £ m (r), increases monotonically with time. 
The precise evolution of £ m (r) with redshift has been studied extensively using both N-body 
simulations (e.g. Jenkins et al. 1998) and analytic methods ( Hamilton et al 1991; Peacock & 
Dodds 1994; Jain, Mo & White 1995). If £ m (r, z) were observable, it would be straightforward 
to use its behaviour to determine Q, A and the power spectrum of linear density fluctuations. 
What one measures in practice, however, is the clustering of galaxies, and the interpretation 
then requires an understanding of how these objects trace the underlying dark matter density 
field. 

If galaxies form at the centre of dark matter halos, considerable insight may be gained by 
using N-body simulations to study the clustering of halos (Brainerd & Villumsen 1994; Mo 
& White 1996; Roukema et al 1997; Jing & Suto 1998; Wechsler et al 1998; Bagla 1998a,b; 
Ma 1998). Mo & White (1996) tested an approximate analytic theory against their numerical 
results, and this theory and its extensions can also be used to analyse the evolution of halo 
clustering with redshift (Matarrese et al 1997; Coles et al 1998). An important conclusion 
from all these studies is that the clustering of halos of galactic mass (~ 1O 12 M ) evolves much 
more slowly than the clustering of the dark matter. This is because at high redshifts, such 
halos correspond to rare peaks in the initial density field, and are thus more strongly clustered 
than the dark matter (Kaiser 1984). Another important conclusion is that more massive halos 
are more strongly clustered than less massive halos. If the luminosity of a galaxy is correlated 
with the mass of its halo, more luminous galaxies ought to be more strongly clustered. A 
detailed comparison with observational data requires a model for the observable properties of 
the galaxies present within halos of given mass at each epoch. 

Precise measurement of the clustering amplitude of galaxies at high redshift has just re- 
cently become feasible. Usually this is done by calculating the angular two-point correlation 
function w(6) as a function of apparent magnitude. In order to assess how clustering has 
evolved, w(6) must be deprojected using Limber's equation under the assumption of some 
specific model for the redshift distribution of the observed galaxies. In future large surveys of 
faint galaxies with photometric and/or spectroscopic redshifts will be available (see, for exam- 
ple, Connolly et al 1995). It will be possible to classify the galaxies in these surveys according 
to absolute magnitude, spectral type, star formation rate and colour, and to investigate how 
clustering evolution depends on these properties. 



1 



At present, most of the data indicate that the clustering amplitude of galaxies decreases 
from z — to z — 1. Different analyses, however, yield very different estimates for the 
strength of this decrease. Le Fevre et al (1996) analyzed the clustering of 591 galaxies with 
/ < 22.5 in the five 10 arcminute fields of the CFRS survey. They find that clustering has 
evolved dramatically, quoting a comoving correlation length at redshift 0.5 of r = 2 h _1 Mpc 
(g = 0.5). More recently, Carlberg et al (1998) presented a preliminary analysis of clustering 
in the CNOC2 field galaxy redshift survey. Their sample is spread over four patches of sky 
with a total area of 1.5 square degrees. They estimate that the comoving correlation length 
of galaxies with M R < -20 evolves as r (z) = (5.15 ± 0.15)(1 + z )-°- 3±0 - 2 h^ 1 Mpc. This is 
much weaker than the evolution found for the CFRS galaxies. The results of Carlberg et al 
(1998) agree reasonably well with those of Postman et al (1998), based on 710 000 galaxies 
with Iab < 24 from an imaging survey of a contiguous 4 square degree region of the sky. The 
latter authors suggest that the small volumes sampled by the CFRS and other early surveys 
resulted in their derived correlation lengths being biased low. 

The clustering of Lyman break galaxies at z ~ 3 has now been measured with surprising 
precision (Steidel et al 1998a; Giavalisco et al 1998; Adelberger et al 1998). The comoving 
correlation length of these objects is comparable to that of galaxies today, implying, as 
expected from models of halo clustering, that Lyman break galaxies are highly biased tracers 
of the dark matter distribution at these redshifts. In addition, Giavalisco et al (1998) find 
that the fainter Lyman break galaxies are less strongly clustered. This accords well with a 
simple model in which the star formation rates in these objects increase with the mass of their 
halos. More detailed theoretical modelling of the observed properties of Lyman break galaxies 
at z — 3, including analysis of their abundances, sizes, luminosities, colours, star formation 
rates and clustering properties, has been carried out by Mo & Fukugita (1996), Baugh et al 
(1998), Governato et al (1998), Somerville, Primack k Faber(1998) and Mo, Mao & White 
(1998) 

In this paper, we combine cosmological N-body simulations and semi-analytic modelling 
of galaxy formation to study the evolution of galaxy clustering as a function of redshift. 
Our methods for incorporating galaxy formation in the simulations are discussed in detail 
in Kauffmann et al (1998, Paper I). Two variants of a cold dark matter (CDM) cosmology 
are analyzed here: a high-density model with Q, — 1, V — 0.2 and H = 50 km s~ Mpc -1 
(rCDM), and a low-density flat model with Q = 0.3, A = 0.7 and H = 70 km s" 1 Mpc" 1 
(ACDM). Paper I was concerned with the global properties of the galaxy distribution at z — 0, 
including B and K-band luminosity functions, the I-band Tully-Fisher relation, galaxy two- 
point correlation functions, colour distributions, star formation rate functions and peculiar 
velocity distributions. Here we focus on clustering evolution in the two models. We study the 
predicted differences in clustering evolution for galaxies of different magnitude, type and star 
formation rate, and we outline how future observational data will clarify the galaxy formation 
process. 
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2 What can be learned from the evolution of halo clus- 
tering? 

Because galaxy formation is complex and involves many poorly-understood physical processes, 
for example, star formation and radiative and hydrodynamical feedback, it is worthwhile to ask 
whether the clustering of dark matter halos can be used to constrain cosmological parameters 
directly. 

In figure 1 we plot the correlation length r as a function of redshift for halos of different 
mass in our two simulations. Here, as in the rest of the paper, all length scales are expressed 
in comoving units. Since the correlation functions in the models are not exact power laws, we 
define r as the radius where £(r) = 1. The smallest halos resolved in the simulations contain 
10 particles and have virial masses ~ 2x 10 n M o . For these objects, r initially decreases 
with redshift, reaches a minimum, and then increases again. The redshift of this minimum is 
different for the two cosmologies: z ~ 0.7 for rCDM and z ~ 1.5 for ACDM. Massive halos do 
not exhibit the same "dip" in correlation length; their remains constant for a while, then 
increases at high redshift. Once again, the redshift at which the evolution becomes strong is 
lower for rCDM than for ACDM. This is simply because structure formation occurs later in 
the rCDM model. 

With 10 metre telescopes, it is now possible to measure the rotation curves of disk galaxies 
at redshift ~ 1 (Vogt et al 1996). It will be many years, however, before such samples are both 
large enough and complete enough for an analysis of the clustering evolution of galaxies as a 
function of their halo mass. In all likelihood, we will have to deal with flux-limited surveys of 
galaxies for some time to come. 

Let us now make the simplifying assumption that each simulated dark matter halo contains 
one observable galaxy, and that the luminosity of the galaxy increases with the mass of its 
halo. The correlation function of a flux-limited sample of galaxies of known abundance at 
redshift z may then be calculated by evaluating £(r) for the mass-limited set of simulated 
halos which has the same abundance. This is illustrated in figure 2, where we plot the 
correlation length of halos versus their number density (in units of h 3 Mpc" 3 ) at a series of 
redshifts. As expected, decreases as the number density increases, because the correlation 
signal becomes dominated by low-mass halos, which are more weakly clustered. Note that the 
differences between the rCDM and the ACDM models are small at all redshifts. Mo, Mao & 
White (1998) show a similar plot for halos at z = 3 for four different CDM cosmologies and 
find that they all give similar results. There thus appears to be a "cosmic conspiracy" that 
makes it impossible to infer information about cosmological parameters from the clustering 
of halos of given abundance. On the other hand, the uniformity seen in figure 2 can be 
used as a test of the entire class of hierarchical models and of the hypothesis that there is a 
one-to-one correspondence between halo mass and galaxy luminosity. As discussed by Steidel 
et al (1998b) and by Mo, Mao & White (1998), this hypothesis works well for the Lyman 
break population at z ~ 3. At low redshifts, the assumption of a one-to-one correspondence 
between halos and galaxies must break down, because the abundance of high mass halos is 
larger, and more and more bright galaxies are grouped together in each such halo. The values 
of r plotted in figure 2 are then lower limits on the true values. 
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Figure 1: The evolution of the co- moving correlation length of halos as a function of redshift 
in the rCDM and ACDM simulations. The solid line is for halos with log (M vir /M & ) in the 
range 11.0 — 11.5, the dotted line for 11.5 — 12, the short-dashed line for 12 — 12.5 and the 
long-dashed line for 12.5 — 13. 
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Figure 2: The co-moving correlation length r of halos is plotted against comoving number 
density at redshifts 0, 0.5, 1, 1.5, 2 and 3. The solid line shows results for the rCDM simulation 
and the dotted line for the ACDM simulation. 
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3 The evolution of galaxy clustering 



In this section, we study the evolution of galaxy clustering in the ACDM and rCDM simu- 
lations. We use the star formation and feedback recipes that resulted in the best fits to the 
observational data at z — 0. As discussed in Paper I, extremely efficient feedback was required 
in the rCDM model in order to obtain a reasonable fit to the correlation function on scales 
below 1 h _1 Mpc and to avoid producing too many galaxies with luminosities below L*. Even 
so, the model failed to fit the observed bright end of the luminosity function and the clustering 
amplitude was too low on large scales. The ACDM model with relatively inefficient feedback 
resulted in a better overall fit to most of the data at z — 0. For simplicity, we do not consider 
dust extinction in the analysis of this paper because it is very uncertain how the empirical 
recipes we adopted in Paper I should be extended to high redshift. This neglect has little 
effect on our ACDM model but means that our rCDM model now substantially underpredicts 
galaxy clustering at z — 0. We concentrate below on the relative evolution of r rather than 
on its absolute value, so this problem does not strongly affect our conclusions. 

In figure 3 results are shown for galaxies with rest-frame B-band magnitudes brighter than 
— 19 + 5 log h in the ACDM simulation. At redshift zero, this corresponds to selecting galaxies 
brighter than ~ L*. The first three panels in the plot show the evolution of £(r) evaluated 
at r = 2, 3 and 8 h -1 Mpc (co-moving units). The fourth panel shows the evolution of the 
co-moving correlation length r . For comparison, the dotted line in each panel shows the 
evolution of the corresponding quantity for the dark matter. Results for rCDM are given in 
figure 4. In each case the redshift extends to the point at which the abundance of L* galaxies 
becomes too low for reliable estimation of the correlation function. 

In the ACDM model, the clustering amplitude decreases from z = to z = 1.5, remains 
approximately constant from z — 1.5 to z — 2.5 and then increases again at higher redshift. 
The dip in clustering amplitude is stronger on small scales: £(r) decreases by a factor of 3 
at 2 h _1 Mpc and by a factor 1.5 at 8 h _1 Mpc. The correlation length tq decreases from 
5.5 h" 1 Mpc at z = to 3.9 h _1 Mpc at z — 1.5. This agrees remarkably well with the 
parametrization of r as a function of z quoted by Carlberg et al (1998). In the rCDM model 
the clustering amplitude remains fixed from z = to z = 1 and then rises steeply at higher 
redshifts. Note, as mentioned above, that the correlation length at z = is low (~ 3 h _1 
Mpc) in this model. 

In figure 5, we plot the evolution of the bias b, defined as the square root of the ratio 
between the galaxy and the dark matter correlation functions: 



The four lines on the plot show the bias as a function of redshift evaluated at r = 2, 3, 5 
and 8 h -1 Mpc. The bias does not depend on r in either model, except at very high redshifts 
where b is somewhat larger on small scales. The bias evolves much more rapidly in the rCDM 
model than in the ACDM model. Galaxies in the ACDM model are unbiased tracers of the 
mass out to z ~ 1. Galaxies in the rCDM model have bias values of about 2 at this redshift. 




(1) 
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Figure 3: Evolution of clustering the the ACDM model. In the first 3 panels, the clustering 
amplitude is plotted against redshift for galaxies with rest-frame B-band magnitude brighter 
than — 19 + 51og/i (solid lines) and for the dark matter (dotted line). Results are shown for 
£(r) evaluated at r = 2, 3 and 8 h -1 Mpc" 1 . In the fourth panel, the comoving correlation 
length r is plotted against redshift both for the galaxies and for the dark matter. 
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Figure 4: As in figure 3, except for the rCDM model. 
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Figure 5: Evolution of the bias for galaxies with rest-frame B-band magnitude brighter than 
—19 + 51og/i in the ACDM and rCDM models. Solid, dotted, short-dashed and long-dashed 
lines show results evaluated on comoving scales of 2, 3, 5 and 8 h _1 Mpc respectively. 

4 Dependence on luminosity, star formation rate and 
morphological type 

In figures 6 and 7, we demonstrate that the clustering evolution depends on the way in which 
galaxies are selected in the simulations. The first panel compares the clustering evolution 
of galaxies selected in the rest-frame B-band with that of galaxies selected in the rest-frame 
I-band. The second panel compares the clustering evolution of galaxies with M(B) < — 19 + 
51og/i with that of galaxies 1.5 magnitudes brighter. The third panel shows what happens 
if galaxies are selected by star formation rate rather than by luminosity. The fourth panel 
shows the clustering of early-type galaxies with stellar mass greater than 3 x 1O 1O M (recall 
from Paper I that these objects form by mergers of two galaxies of similar mass and have 
M(B) bulge - M(B) total < 1 mag). 

In the ACDM model, we find that the strength of the "dip" in clustering between z = 
and z = 1.5 is sensitive to sample selection. The dip is weaker for more luminous galaxies and 
for galaxies selected in the I-band, but stronger for galaxies selected by star formation rate. 
The clustering amplitude of early-type galaxies is stronger than that of the population as a 
whole and evolves very little with redshift. In the rCDM model, clustering evolution is less 
sensitive to sample selection. The clustering always remains fixed out to z ~ 1 and then rises 
steeply at higher redshifts. Early-type galaxies are also very strongly clustered in this model, 
particularly at high redshifts. 
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Figure 6: The dependence of clustering evolution on sample selection in the ACDM model. 
The comoving correlation length r is plotted as a function of redshift for: a) L* galaxies 
selected in the rest-frame I-band (solid) and galaxies selected in the rest-frame B-band 
(dotted); b) very bright galaxies (solid) and galaxies (dotted); c) galaxies with star forma- 
tion rate greater than 3 M & yr" 1 (solid) and 1 M yr _1 (dotted); d) early-type galaxies with 
stellar masses greater than 3 x 1O 1O M . 



10 




11 



4.1 What can be learned from these dependences? 

We now explain why the evolution of clustering depends on sample selection and what can be 
learned about galaxy formation by studying the observed evolution as a function of morpho- 
logical type, luminosity, colour and star formation rate. 

For a given cosmology, the clustering amplitude predicted for a sample of galaxies depends 
on the masses of the dark matter halos they inhabit. The evolution of clustering depends 
on how the mass distribution of these halos changes with redshift. Additional relevant and 
observationally accessible information comes from the variation with redshift of the abundance 
of galaxies in the sample. 

As an example, let us suppose that galaxies with fixed star formation rate are found 
in smaller halos at high redshift than at the present day. We would expect the clustering 
amplitude of a SFR-selected sample to show a stronger dip than a sample of galaxies that 
tracked halos of the same mass at all redshifts. We would also expect the abundance of 
galaxies in a SFR-selected sample to increase more strongly with redshift, because there are 
many more small halos than large ones. 

As a second example, let us suppose that early type galaxies are found primarily in massive 
halos at all redshifts. As seen in figure 1, these galaxies should not exhibit any dip in clustering 
and their abundances should decrease strongly at high redshifts because massive halos are rare 
objects at early times. 

These points are illustrated in detail in figure 8, where we plot the evolution of the median 
halo mass and the comoving number density of galaxies in samples selected in different ways 
from the ACDM simulation. The top 3 panels show results for galaxies selected according 
to rest-frame B- magnitude, rest-frame I- magnitude and star formation rate. At z — 0, all 
three galaxy samples have the same abundance and occur in halos of roughly the same mass. 
Galaxies selected according to star formation rate move to smaller halos at higher redshift. 
This effect is simply a result of the parametrization of star formation in our models. Following 
Kennicutt (1998), we have adopted a star formation law of the form M* = aM co id/td yn , where 
M co id is the mass of cold gas in the galaxy and td yn is the dynamical time of the galaxy. Since 
tdyn decreases at higher redshifts, the star formation rates are higher in halos of the same 
cold gas content. Galaxies selected in the B-band exhibit a weaker trend towards low- mass 
halos. In the case of the I-band selection, galaxies trace halos of roughly the same mass at 
all redshifts below 2. This is because the I-band magnitude of a galaxy is a measure of its 
total stellar mass, rather than its instantaneous star formation rate. We thus conclude that 
SFR-selected samples show the strongest dip in clustering in figure 6 because this selection 
procedure favours galaxies in lower mass halos at high redshift. Note that galaxies in the 
SFR-selected samples also exhibit the strongest increase in abundance from z — to z — 1.5. 

The bottom two panels in figure 8 show that very bright galaxies and early-type galaxies 
in the simulation are found in halos with masses ~ 10 13 M Q . As seen in figure 1, the clustering 
of these objects evolves very little from z — to z — 1. 

In figure 9, we show the evolution of galaxy abundances and median halo masses for 
samples selected in the same way from the rCDM simulation. The results are qualitatively 
similar to those found for ACDM. galaxies occur in halos of roughly the same masses 
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Figure 8: The evolution of the comoving number density (left column) and the median halo 
mass (right column) of galaxies selected from the ACDM simulation. Error bars indicate the 
upper and lower quartiles of the halo mass distributions. Results are shown for the selection 
criteria described in the caption to figure 6. 
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(~ 1O 12 M ) in both models. The reason why no dip is seen in the rCDM model is because 
halos of these masses are more strongly biased at z — 1 than in the ACDM model and the 
decrease in halo mass with redshift for the SFR-selected sample is less pronounced. Note also 
that the redshift at which the abundance curves peak is higher for ACDM than for rCDM. In 
the ACDM simulation, the abundance of early-type galaxies only decreases substantially at 
redshifts greater than 1.5, whereas in the rCDM simulation, the abundance of ellipticals has 
already declined by a factor of 3 by z — 1. 

5 Evolution of the slope of £(r) 

Figure 10 shows the evolution of the slope 7 of the two-point correlation function for galaxies 
with rest-frame B-magnitudes brighter than —19 + 51og/i in the ACDM and rCDM models. 
We have fit a power-law to £(r) over three different ranges in scale: r = 1 — 5 h _1 Mpc, 
r = 5 - 10 h" 1 Mpc and r = 1 - 10 h" 1 Mpc. 

In the ACDM model, the evolution of the slope is stronger on small scales. Over the range 
1—5 h _1 Mpc, 7 evolves from —2 at z = to —1.5 at z — 1. On large scales, 7 remains 
approximately constant. Over the range 1 — 10 h _1 Mpc, 7 evolves from —1.85 at z = to 
—1.6 at z — 1 and then remains constant. These results appear to be in qualitative agreement 
with the observations (Postman et al 1998). These au thors find no dependence of 7 on 
magnitude for the bright (/ < 21) galaxies in their survey. At fainter magnitudes 7 flattens, 
reaching a value of —1.6 at 1=22.5. They also find that the flattening is stronger on smaller 
angular scales. Neuschaeffer & Windhorst (1995) find similar results from an independent 
survey carried out at a different wavelength. In the rCDM model, there is very little change 
in the slope with redshift on any scale. 

6 Evolution of pairwise peculiar velocities 

Estimates of dynamical quantities such as Q or cluster M/L ratios from galaxy data require 
knowledge not only of the spatial bias in the galaxy distribution, but also of any possible bias 
in the kinematics of the galaxies relative to those of the dark matter. In Paper I we explored 
this " velocity bias" for our z = models using pairwise velocity statistics, and in Paper III 
(Diaferio et al 1998) we will do the same using group and cluster velocity dispersions. Here we 
briefly explore the predicted evolution of velocity bias using pairwise statistics. The thin solid 
lines in figure 11 show the redshift evolution of the pairwise peculiar velocity dispersion a 12 
evaluated at relative separations r = 0.5, 1 and 2 h _1 Mpc (co-moving units) for galaxies with 
rest-frame B-band magnitudes less than —19 + 51og/i in the ACDM and rCDM simulations. 
The thick solid lines show the evolution of a 12 for the dark matter. In order to compare the 
relative change in o\ 2 as a function of redshift in the two models, we scale the results by 
dividing by the value of a{ 2 at z=0. As shown in figure 13 of Paper I, af^ — 800 km s _1 
(r = 1 IT 1 Mpc) in both the ACDM and rCDM models at the present day. 

The pairwise peculiar velocities of the galaxies follow those of the dark matter quite closely 
in both models. The galaxy velocities are 10-40% lower than those of the dark matter at 
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Figure 9: As in figure 8, except for the rCDM simulation. 
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Figure 10: The evolution of the slope of the correlation function of galaxies with M(B) < 
— 19 + 5 log h in the ACDM and rCDM simulations. The solid line is the result of a fit to £(r) 
over scales between 1 and 5 h -1 Mpc, the dotted line is for scales between 5 and 10 h -1 Mpc 
and the dashed line for scales between 1 and 10 h" 1 Mpc. 
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Figure 11: Redshift evolution of the pairwise velocity dispersion a 12 at relative comoving 
separations of 0.5, 1 and 2 h -1 Mpc for galaxies with rest-frame B-magnitudes brighter than 
— 19 + 51og/i (thin lines) and for dark matter (thick lines) in the ACDM and rCDM simula- 
tions. The results are scaled by dividing by the value of a{f at z — 0. 

z = 0. (Note that the "antibias" in galaxy peculiar velocities at z = is stronger than that 
shown in figure 13 of Paper I, because the models presented in this paper do not include dust 
extinction. Dust reduces the contribution of star-forming field galaxies in a B-selected sample, 
but has little effect on early-type galaxies in rich groups and clusters. Models that include 
dust extinction thus give values of afjf that are 10-25 % larger ). The difference between the 
galaxy and dark matter peculiar velocities decreases at higher redshift. In contrast to the 
spatial distributions, galaxy peculiar velocities in our models are never very strongly biased. 
In the rCDM model, there is nearly a factor 2 decrease in a{f from z=0 to z=0.5. In the 
ACDM models, ofjf remains roughly constant out to z = 0.5, before decreasing at higher 
redshift. 



7 Discussion and Conclusions 

In a hierarchical Universe, the evolution of galaxy clustering depends on the following: 
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1. Cosmological parameters, such as Q, A and a$, because these determine the rate at 
which structure grows and the epoch at which halos of given mass change from being 
rare objects, and thus biased tracers of the dark matter distribution, to being "typical" 
objects with clustering properties similar to that of the mass. 

2. The relationship between the mass of a dark matter halo and the properties of the 
galaxies that form within it. Note that this relationship depends only on halo mass and 
is independent of the environment in which the halo finds itself (Lemson & Kauffmann 
1998). 

3. The evolution of the galaxy population with redshift (and hence the evolution of the 
relationship between galaxy properties and halo mass). 

In this paper, we illustrate how differences in clustering evolution between galaxies of 
differing luminosity, colour, morphological type and star formation rate may help constrain 
galaxy formation models, and perhaps even cosmological parameters. 

One interesting diagnostic that we highlight is the "dip" in correlation amplitude observed 
between z = and z — 1. We show that this dip occurs naturally in a ACDM model, 
where structure forms early and halos with masses in the range 10 11 — 1O 12 M , which contain 
galaxies of intermediate luminosities, are unbiased tracers of the mass over this redshift range. 
In the rCDM model, bias evolves rapidly, and the clustering amplitude of L* galaxies remains 
constant from z — to z — 1. Although it might be possible to "force" a dip in clustering in 
the rCDM model by requiring that galaxies form in less massive halos, it would be difficult 
to come up with a physically-motivated scheme for doing this that would not simultaneously 
produce too many bright galaxies. 

We also show that the strength of the dip in the ACDM model is sensitive to sample selec- 
tion. If galaxies are selected according to star formation rate rather than B-band luminosity, 
objects in low mass halos contribute more to the clustering signal at high redshifts and the 
dip is stronger. If galaxies are selected in the red rather than the blue, the dip is reduced. 
Very luminous galaxies and massive early-type galaxies exhibit no dip in clustering between 
z = and z — 1 because they occur in high mass halos (10 13 — 1O 14 M ) that are already 
biased at z = and become substantially more biased at high redshift. 

The predictions presented in this paper should be viewed as illustrative rather than quan- 
titative. As discussed in Paper I, the precise relation between the mass of a halo and the 
properties of the galaxies that form within it depends strongly on the adopted recipes for 
star formation and feedback; these are very uncertain. The exciting prospect is that future 
observations of galaxy clustering at high redshift will place strong empirical constraints on 
these processes. 
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